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INTRODUCTION 


A . Motivation 

In the discipline of microwave radiometry it is often important to determine the 
complete emission characteristic of a flat scene. The emission characteristic is commonly 
estimated from a set of antenna temperature measurements conducted at a number 
of observation angles from a remote observation point. It is well known that good 
estimates of the brightness temperature distribution require that the measurements 
be compensated for (1) antenna pattern weighting of the distributions, and (2) direct 
and reflected atmospheric emission contributions. Compensation for antenna pattern 
weighting yields the apparent temperature distribution about the point of observation. 
Extraction of the direct and reflected contributions from the apparent temperature 
distribution then yields the brightness temperature characteristic of the scene. 

Corrections of either type represent rather formidable tasks both in the implementation 
of the experiment and in the reduction of the data. The precision inherent in the 
radiometer can, however, only be realized in estimating the temperature distribution 
when thorough and efficient methods for handling both types of corrections are 
employed. 

Various methods to estimate the apparent temperature have been used. In 
some instances it was common practice to use the antenna temperature as an estimate 
of the apparent temperature (and sometimes mistakenly as the brightness temperature). 
However, Peake et al .[ 1 ] , using a recursive inversion technique, have vividly 
illustrated the difference between the antenna temperature and the apparent temperature. 
Fourier transform techniques have been reported [2,3,4], largely by members of the 
radio astronomy community. In these methods repeated and sometimes lengthy numeri- 
cal double integrations are required to estimate the apparent temperature for each 
scene. They are, as a consequence, inefficient for handling large volumes of data 
as often collected by those engaged in remote sensing programs. 

An efficient as well as comprehensive method for estimating the apparent 
temperature distribution has been derived for isotropic scenes. The ahtenna weights 
on the vertically and horizontally polarized emission components are formed by a 
linear mathematical model which when invertible yields the apparent temperature 
distribution. It is shown that radiometric observations at both polarizations are 
required at a number of angles between nadir and zenith to recover the apparent temp- 
erature distributions. 


B. Background 

The antenna temperature inversion problem is mathematically akin to a host 
of physical problems in which a function f(x) is to be inferred from measurements 
g(y) taken by a linear system where the measurements are related to f(x) by a 
Fredholm equation of the first kind 

5<3 > = f b K < *, 3 ) f <*> d* ( 1 ) 

J ft 

In this particular application f(x) is to be regarded as the apparent temperature 
distribution, the kernel K(x,y) as the antenna power pattern and g(y) as the antenna 
temperature . 

Various inversion techniques for (1) have been proposed. When main beam 
smoothing is the dominant effect, a "desmoothing" technique as suggested by 
Bracewell [ 5 ] can be employed. In the general inversion problem, more than the 
main beam will play an important role, however. When the functional form of f(x) 
is known on an a priori basis, except for a set of unknown parameters, very often the 
optimum member of a family can be determined by comparing the measurements in 
the least squares sense with an algebraic representation of (1). This type of approach 
is often suggested for sky temperature measurements [6]. However, for terrestial 
and ocean scenes the apparent temperature distribution is not known on an a priori 
basis. Furthermore, any technique suggested for sky temperature measurements would 
not be appropriate for the general case since emissions from the atmosphere are 
usually considered independent of polarization. An iterative technique similar to 
to Picard's process of successive approximations [7] has been reported [1,8]. Sometimes 
called the boot-strap technique, this method uses the measured g(x) as the first 
estimate of f(x). When the measurement is subject to noise, convergence of the 
iterative technique is not assured in practice [8], [9]. A Fourier inversion technique 
based on interferometric measurements is employed by the radio astronomer. The 
technique, however, is restricted to sources of small angular extent [3], [10], Here 
a method is sought for expansive flat scenes such as the ocean. 
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C . Approach 

A matrix formulation of the integral equation relating the apparent temperature 
to the antenna temperature is employed here. Pictorial ly, the antenna is located at 
the center of a spherical coordinate system and is observing a cylindrical ly symmetric 
temperature distribution* on a sphere. The cylindrical symmetry applies to the upper 
and lower half spaces separately. The temperature distribution, however, is approxi- 
mated by spherical bands at constant temperature as illustrated in Figure 1. When 
the antenna "points" at the midpoint of one of these bands, the antenna temperature 
is considered a linear sum of the antenna weights on each band. The antenna weights 
for observations on all bands form the modeling matrix. This approach must therefore 
assume that the apparent temperature distributions are smooth functions so that they 
can be approximated by a constant over each band.** The antenna main beam is 
considered sufficiently narrow to spatially resolve temperature differences between 
adjacent bands. 

Particular attention is given to generalized descriptions of the antenna patterns. 
Here, it is assumed that the patterns were measured from an azimuth-over-elevation 
mount so as to preserve the notion of polarizations as defined in a spherical coordinate 
system. Each polarized pattern is considered to be accompanied by a cross-polarized 
pattern of generally smaller intensity. The cross polarized component is assumed to 
be coupled in phase with the polarized component. The antenna polarizations (some- 
what surprisingly) do not in general coincide with the surface polarizations.*** 

The formulation as a consequence distinguishes the two carefully and therefore 
accounts for cross- coup ling between antenna and surface polarizations. 

In a matrix formulation of this type the user must be keenly aware that under 
some circumstances the linear model can become ill-conditioned [11], particularly 
when the model is developed at the limit of the antenna resolution. In some 
circumstances an ill-conditioned model can be salvaged by amending the formulation 
with numerical filtering [12], 


* Not a necessary assumption, although once discarded the formulation becomes 
unwieldy as well as the measurements. 

** Experience has shown that several bands can violate this assumption without 
adversely affecting the inversion estimates on the other bands. 

.his pecuSarity was rnaue ciear Dy vv . n. reaxe or xumo orare university in a 
private communication. 
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THEORETICAL DEVELOPMENT 


A. Description of the Antenna Pattern 

A dual polarized antenna is located at the origin of a primed coordinate system 
(see Figure 2). The boresight is directed along the positive x' axis. When the v~ 
port is energized and the h port is terminated, the vertically polarized pattern is 
generated. The E plane coincides with x'z' plane. Similarly, when the h port is 
energized and the v port terminated, the horizontally polarized pattern is generated. 
The x'y' plane coincides with its E plane. Each polarized pattern is assumed to be 
accompanied by a cross-polarized pattern. In general, the polarized field is coupled 
in phase to the cross - polarized field when the antenna is regarded as a transmitter. The 
resultant field is consequently elliptically polarized. By the reciprocity theorem it 
is understood that the reception property has a similar polarization character. The 
phase coupling, however, must be regarded as the complex conjugate of that when 
the antenna is viewed as a transmitter. For the polarization discriminant antenna, 
the effect is weak since the cross polarized pattern is assumed to have considerably 
smaller power content. However, a general viewpoint will be retained since the 
radiometer antenna receives partially polarized fields [13], The orientations of the 
antenna polarizations are defined in terms of spherical unit vectors as illustrated in 
Figure 2. 

In view of the elliptically polarized character of the antenna, the 

reception property will be described in terms of a vector complex effective 

th 

height. For the p (p = v or h) port the effective height is denoted by 


a 


P 


/4tt 



<e>'> 


L? 


( 2 ) 
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where 


A = 

G (eT) = 


mean operating wavelength 
elliptically polarized power gain 


f P 


( 3 ) 


R may be identified with the reception factor defined by Kales [14]. In expression 

P 

(3) Y and 5 are complex functions of 9" and <t>' and denote the relative sensitivity 

P P 2 

to a polarized field and a cross polarized field. The components of | Rp| are 

related to the typically measured normalized patterns g ] and g « by the following 

P P^ 

expressions 


3 P i 
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A 'Tr 


( 4 ) 
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where 

9 vl <3h2 ) = 


vertically (horizontally) polarized pattern 


« /_ \ - 
a v2 vy hl' 


cross polarized pattern associated 

whh g vl (g h2 ) 


7 


) d-a 


(6) 


•a-t * JJ ( 3f + 9p2 

-Hr 

H = radiation efficiency 

P 

In the following development it will be convenient to consider the radiation efficiency 
factor as unity and to incorporate its effect in the radiometer transfer function. The 
antenna is therefore idealized. 

The relative phase between the polarized and cross polarized fields must be 

measured to establish the phasor characters of Y and 6 . To understand the nature 

P P 

of P it is helpful to write it in the coordinate frame coinciding with the principal 
axes of the polarization ellipse. Let ij and i ^ denote orthogonal unit vectors coin- 
ciding with the principal axes for a particular look direction. Then 


= cos % i., + j sin % Lz 


(7) 


where | tan X | is the axial ratio. Now let £ denote the angle of rotation between 
the principal axes and the antenna polarization vectors as illustrated in Figure 3. 
Then 


f = ( cos X cos Si - j sLn X. sin *■ ) 

« 

+ ( cos X stia ^ -*■ j sen X cos ^ (8) 


The terms within the brackets are identified as Y 
particular look direction. 


P 


and <5 , 

P 


respectively, for this 


B . Geometry of the Measurement Coordinate System 

The above described antenna is now positioned at a height H above a "flat" 
scene whose emission characteristic is to be measured. The antenna is considered to 
be located in an absorbing/emitting atmosphere. The geometry of the scene is 
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described in terms of an unprimed coordinate system whose z axis points downward 
(see Figure 4). The antenna scans in the xz plane from nadir to zenith. In this 
plane the look direction 0 O is defined as the angle between the boresight axis x' 
and the z axis. The angle 6, also describes the rotation about an axis common to 
the antenna and scene coordinate frames. 

C . Development of the Inversion Model 

The antenna observes an angular distribution of non-coherent sources represented 
by the electric field intensity 


; 8 + E .1* 
— , — 


(9) 


where Eg and are complex functions of 0 and t. Since the bandwidth of receiver 
B is considered much smaller than the frequency of observation,both Eg and E^ may 
be considered quasi - monochromatic in character [ 13 ]. The emission components 
of the surface below or, for that matter, of a sphere surrounding the observation 
point are described in the unprimed coordinate system. The emitted polarizations 
are so chosen to agree with the standard definitions for the polarizations of a planar 
surface (see Figure 4). For an observation in the 6 0 direction the. average power 
available at the p*^ 1 port of the antenna is therefore given by 


p f - ^^7 JJj i e •«,'*> 

p * v or h 


•*lr 


( 10 ) 


where 
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FIGURE 4. GEOMETRY OF THE OBSERVATION COORDINATE SYSTEM 

RELATIVE TO THAT OF THE ANTENNA COORDINATE SYSTEM 
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and 


intrinsic impedance 


n = 

o 


The integrand of (10) when expanded, contains a number of terms involving 

★ * 

auto- correlations and cross-correlations of the type<Eg Eg> and<E^ / respectively 


~<P G' 


For a flat scene it is easily shown that the cross- correlation terms vanish when 
the emissions are formulated in terms of the vertically ig and horizontally ^po larized 
states. On the otherhand, the auto - corre lotion terms are related to the apparent 
temperatures Tq and by the expressions [15] 




A Tr l< T. B 


(12) 


and 


< | E*J> _ ±jrk T» B 


r (« 


(13) 


Here k is Boltzmann's constant. Tg is the vertically polarized contribution; 
whereas is the horizontally polarized contribution. 

The vector inner products generated in the integrand may be reduced to 
trigonometric expressions. It is noted that the orthogonal pairs (i ,i . ) and 
(tg 1 ,^ ')are mutually orthogonal to a line of sight vector The orthogonal 
pairs are therefore related by a simple rotation as illustrated in Figure 5. It follows 
that 
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© 

II 

?*• 

* if 

- cos 


-c <±* * - 

* if 

= sin y 

(14) 


When the received power P ( 0 O ) is translated into an equivalent antenna 

P 

temperature. Equation (10) in view of the above observations and equations (4) and 
(5),can be written as 
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FIGURE 5. RELATIONSHIP BETWEEN ANTENNA AND 
SURFACE POLARIZATIONS 
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where 


t, <e-> * -X7 If { T ° [ 3 ,< * 3 f' sin ‘+ J 

, t < 15) 

+ T f ( 3 pl sin f + 3 ps cos f ) 

+ C T, - T> ] J 3f , 3r ] cos p t Sin 2 t } j« 
* 

^ U.c** p *>J .... 


Suppose the temperature distributions are sampled at both antenna polarizations 
over the interval 049 o < 7 r. Denote the observation angles as (X0j<e < . . . 

Under the assumption that the temperature distributions can be approximated 
by spherical bands at constant temperatures, the antenna temperatures at 6^ may be 
approximated by 


T » l V 

* r 

[% T . ( v + \ 

T * <v 

] 

07) 

T k <e k) 

■n 

• c 

i »« 

Is T » < e j > ♦ Jkj 

T »<v 

) 

(18) 


where 


\ ‘a ; // f 3». cos ** ’ 9„ f /iX cos /*« ] dn 

k '‘j ■ a: jy f s vp sin '' t ' + 9„ c ° s, t - y 9 v, 3 , £ 3111 cos / 3 v ] <*■"• 

(19) 

%■ “ AT //[ 3k. “ s *t + 3k, sU V/3k,3k, sinifcos^J jfl 


d kj “‘X7 / /f3 l( * in * t * 3 k , - J 3 k, 3 k j ^^cosp k ) ia 
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and 


il. = 


r e i-, + ®i ^ + 0 j-. 

{<«,*> • J 2 " 3 " * e * - " 2 — , 

0 < <f> < Z Tr f | < j < n j 

{ (e,f) : 0 < 0 < t>,/2 , 0 < <f> < 2 tt ]■ 




- {ce,*: 


7T+ 


<e < -Jr , 0 < ^ 2Tr } 


The integrations are to be specifically performed in the unprimed coordinate system 
as defined by ft.. As will be seen, the integrands are functions of 9^ and thereby 
justify the kj subscripts for the constants a, b, c and d. The density of measure- 
ments must be sufficient to justify the approximation that the apparent temperatures, 
Tq and T^, are essentially constant on each spherical band ft .. 

The above system of equations (11) and (12) can be written in a compact form 
with partitioned matrices 
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and (A ) 
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<x, 


(B VS 


<C V c ‘j 


<»\j ■ J 


(23) 


The matrix [A,B] contains the antenna weights for the vertically polarized observations 
and [C, D] contains the antenna weights for the horizontally polarized observations. 
The B matrix denotes the horizontally polarized contributions to the vertically polarized 
observations, whereas the C matrix denotes the vertically polarized contributions to 
the horizontally polarized observations. When the modelling matrix is invertible, 
the apparent temperature distribution is approximated by 
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(24) 
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-1 
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C D 
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For symmetric and identical v and h antenna patterns the modeling matrix . (21) 
must be simplified when observations are conducted at nadir and zenith. The first 
and last rows of [A B] are, in this case, identical to the first and last rows of [C D] . 

The isotropic assumption permits removal of the redundant rows, say those inJC D] . 

The corresponding antenna weights located in the_first and last columns of /however, 
must be added to the first and last columns of £ j /respectively. The resulting model 
is then expressed as 


T v 


A 

C' 


B 

9 

D 


T* 

lVj 


(25) 


where primes reflect the removal of rows and columns and the summation of the 
indicated columns. 

It is of interest to note that for a non-precipitating atmosphere the emissions 
coming totally from the atmosphere are polarization invariant. As a consequence 
the latter elements in Tq , say Tg (0^ )/ T e ^k+i) • • • ' T 0 (Q^' are redundant 
with the respective latter elements in Tj, , say T^ (0^), T^ (9^+]) • • • / (9 n ). In 

this circumstance, the modelling matrix can be altered in either of two ways. If 
observations above the horizons are conducted only at one polarization, the modelling 
matrix may be collapsed by removing the rows associated with the deleted observations. 
The associative columns must be removed and added to the columns corresponding 
to the polarized observations which have been retained. The result is a square matrix 
whose inverse will estimate the apparent temperature distributions. On the otherhand, 
if observations at both polarizations are retained above the horizon, the original 
modelling matrix may be solved subject to the constraint that Tg (8,) = ^ (0.), j=k, 
k+1, . . ., n. In this case, all the rows are retained; however, columns k through 
n in are removed and added to the corresponding columns of gj • The result- 

ing rectangular system may be solved by the least squares method. It is anticipated that 
the latter technique will yield better temperature estimates above the horizon than 
the initial one. 
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Under the isotropic assumption either one of the above two procedures (applied 
to the observations above the horizon) may also be applied to the nadir observations. 
The latter procedure will also work for symmetric patterns. 


D . Computation of the Elements in the Model 

From a practical viewpoint the elements in the radiometer inversion model 
(21) must be evaluated numerically. Since the integrations are to be performed 
over spherical bands the arguments of the integrand functions, as a matter of 
convenience, must be expressed in terms of the unprimed coordinates. The transformation 
between the coordinate systems is clearly (see Figure 4) given by 
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Now from (14) it is recalled that 
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(32) 
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has been used. 

Relationships (29) and (31 ) resolve the ¥ and 4> ' dependence on 9 and . 
The 0' dependence on 0 and <t> is established by 


or 


cos e' = 


-t 

— r 
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— i 


(33) 


cos e' = cos e si. - *, e, - si.n & cos e 0 cos f (34) 


SIMULATION RESULTS 

A. Introduction 

To evaluate the radiometer inversion model posed above, computer simulation 
studies were conducted. With an assumed set of antenna patterns and assumed 
apparent temperature distributions, antenna temperatures were accurately computed 
at ten degree intervals. These antenna temperatures were employed in a matrix 
inversion model based on observations at 10 degree intervals. By means of the inverse, 
the apparent temperature distributions were recovered and compared with the original 
temperature distributions. The implementation and results of these studies are 
described below. 

B. Computation of Elements in the Model 

Functional representations of the AAFE RADSCAT antenna* patterns were 
employed to evaluate the weighting factors (elements) in the radiometer inversion 
model (17) - (20). The h plane pattern at 13.9 GHz is shown in Figures 6. The 
polarized main beam, the first few side - lobes (partially embedded in main beam), 
the cross-polarized beams, and the spill-over side lobes, were all approximated with 
Gaussian functions. The band of side lobes around the main beam was approximated 
with a cylindrical cone as shown by the dashed lines in Figure 6. Since the back lobes 
were insignificant, they were ignored in these studies. In addition to the polarized 
pattern, a representation for the cross polarized pattern (not shown) was also employed. 

TC 

The AAFE RADSCAT antenna is a rear fed parabolic antenna developed under 
NASA Contract NAS 1-10161 by General Electric Space Division, Philadelphia, Pa. 
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The cross polarized pattern was primarily significant in a vicinity about boresight. 

The relative phase 3 p between polarizations was not available during these 
simulation studies. Fortunately for this antenna, the factors in the inversion model 
involving 3 may be regarded as negligible since the cross polarized pattern was 
localized about the boresight site axis and its maximum intensity was more than 20 dB 
beneath the polarized pattern. This implies Vg pl g p2 sin 24* is insignificant in 
comparison to other terms in the vicinity of the boresight axis. (However, for another 
antenna this factor may be significant). The pattern representations should be considered 
expedient ones which are adequate to evaluate the inversion model rather than 
accurately specifying the actual patterns. The major effects have been preserved 
however . 

The pattern representations together with the polarization mis-match factors 

2 2 

cos y and sin were employed to evaluate numerically the elements a^., b^., 
c^., and d^. in the model. Numerical integrations over the spherical bands ft. 
were performed with Gaussian - Legendre quadrature techniques [16]. The numerical 
integrations roughly require the same amount of computer time regardless of the number 
of observation points { 9 ^: k= l,2,...,n} since it is governed primarily by the 
antenna pattern sampling requirement of the quadrature method. Numerical integration 
in two or more dimensions typically consumes large amounts of computer time and, in 
this instance, since a transformation is required between the domain of integration 
and the pattern coordinates, the time is substantially larger. 

The relative antenna weighting factors on 5° bands for a measurement at an 
observation angle ( 0 ) of 40° on the v port are illustrated in the graphs of Figure 
7. The solid line, plotted from the A matrix, represents the antenna sensitivity to 
the vertically polarized emissions; whereas the dashed line, plotted from the B matrix, 
represents the antenna sensitivity to the horizontally polarized emissions. The 
graphs demonstrate that the antenna temperature will be composed of contributions 
from both surface polarizations virtually over the whole space with the primary part 
coming from the main beam. Although no band off the main beam contributes 
significantly when considered alone, the total contribution from all these bands can 
be significant (even for this reasonably efficient antenna). The largest weight 
on the horizontally polarized emission (nearly 1%) is largely attributable to the 
cross polarized pattern. On the otherhand, the remaining horizontally polarized 
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weights are entirely a result of cross coupling (g y -j sin V ) between the polarized 
I pattern and the "cross" polarized surface emission since the cross polarized pattern 

vanished outside the immediate vicinity of the main beam. From the graphs of 
Figure 7 it is observed that this cross coupling effect is comparable with the like 
coupling effect for contributions entering off the main beam. 

A 38 x 38 modeling matrix representing 19 sample points at 10° intervals 
i and both polarization were prepared. However, the inversions for the apparent 

| temperatures were actually performed with a 38 x 28 matrix which resulted from 

v considering the apparent temperature polarization invariant above the horizon and 

at nadir. 

C . Computed Antenna Temperature Characteristics 

The antenna temperatures at 10 degree intervals were accurately computed 
for temperature distributions typifying the emission characteristic over the sea as 
observed at 7 kilometers for two sky conditions. The distributions for the two sky 
conditions are shown in the graphs of Figures 8 and 9 and are based on the efforts 
of Wu [17] and Wu and Fung [18]. The differences between the antenna temperatures 
and actual polarized temperature distributions are shown in Figures 10 and 11 . It 
is readily apparent that the antenna temperature is a crude measure of the actual 
i temperature distribution. The deviation from the apparent temperature is typically 

| largest at 90 degrees. Comparative studies of the polarized antenna temperatures 

would be sensitive to the different error characteristics as demonstrated by the graphs. 

D . Estimating the Apparent Temperature 

The above simulated antenna temperatures were used in the rectangular inversion 
model described above. The least squares solutions yielded estimates of the apparent 
temperatures for both sky conditions. The differences between the computed and actual 
’ temperatures are shown in Figures 12 through 15 for the two sky conditions and two 

polarizations. From these graphs it is apparent that better than an order of magnitude 
of improvement over that provided by the antenna temperatures can be realized for 
most look angles. The error in the temperature estimates for 0 o >14O° are not shown 
in the graphs since the errors were negligible. In contrast to the antenna temperature 
error characteristics, the computed apparent temperature errors are largely independent 
of polarizations and the actual distributions. It is anticipated that the errors near 
the horizon could be reduced by sampling there more often. 
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FIGURE 7. MODEL WEIGHTING FACTORS AS A FUNCTION OF POLAR ANGLE 

FOR A VERTICALLY POLARIZED OBSERVATION AT 40 DEGREES 
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FIGURE 9. TEMPERATURE DISTRIBUTIONS SIMULATING HEAVY OVERCAST 

CONDITIONS OVER THE OCEAN AS OBSERVED AT AN ALTITUDE 
OF 7 KILOMETERS 
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FIGURE 10. COMPARISO 
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FIGURE 11. COMPARISON OF THE ANTENNA TEMPERATURE WITH THE RESPECTIVE 
POLARIZED TEMPERATURE DISTRIBUTION FOR THE OVERCAST CASE. 
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FIGURE 12. COMPARISON OF COMPUTED AND ACTUAL APPARENT 
TEMPERATURE FOR CLEAR SKY CASE. 
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FIGURE 13. COMPARISON OF COMPUTED AND ACTUAL APPARENT 
TEMPERATURE FOR CLEAR SKY CASE. 
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FIGURE H. COMPARISON OF COMPUTED AND ACTUAL APPARENT TEMPERATURES 
FOR OVERCAST SKY CONDITION. 
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FIGURE 15. COMPARISON OF COMPUTED AND ACTUAL APPARENT 
TEMPERATURES FOR OVERCAST SKY CONDITION. 


CONCLUSIONS AND COMMENTS 


An integral inversion method based on a matrix approximation of the integral 
has been successfully applied to yield good estimates of the apparent temperature 
distributions from simulated antenna temperatures. These results have shown that 
cross coupling between antenna polarizations and surface polarizations contributes 
significantly to the antenna temperature. Any inversion method for scenes having a 
polarized character must, as a consequence, account for contributions from both 
emitted polarizations. To adequate ly discern the polarization effects in the inversion 
model, comprehensive knowledge of the polarized and cross _ polarized antenna 
patterns are also required to compute the entries in the observation matrix. In view 
of the cross _ coupling and the existence of cross polarized patterns, it is necessary 
that measurements be conducted at both polarizations to realize a thorough inversion. 

The sensitivity of the above model to errors in the observation vector 
or in the modeling matrix were not considered above. The reader may astute I 
question the behavior of the solution technique in the presence of errors. Some 
preliminary results have indicated that the errors in the estimate follow the error 
character inserted in the observation vector. The dominant diagonal of the modeling 
matrix would lead one to expect this result. Numerical filtering [12] reduced 
the error sensitivity slightly. Sensitivities in the modeling matrix have not yet 
been examined . 

Although the preparation of the inversion model represents a sizeable computer 
effort (in addition to the measurement of the patterns), once the model has been prepared, 
it may be rapidly and repeatedly applied to observations from many scenes which have 
been sampled at the same intervals. The method, in this respect, is efficient for large 
volumes of data, when compared to iterative or Fourier transform methods where 
integrations are required repeatedly for each scene. 

These simulation studies and studies performed by Peake et al. [1] have 
demonstrated that antenna temperatures serve as a crude and sometimes poor measure 
of the actual temperature distribution. The problem is particularly acute when side, 
spill - over, and back lobes make significant contributions to the antenna temperature. 

Thus precision radiometry, in the sense of recovering the apparent temperature 
distributions, involves more than precision at the antenna terminals. Unless highly 
efficient, narrow~beam antennas having low cross - polarized patterns are employed, 
compensation for antenna pattern weighting effects is clearly a prerequisite for 
exacting studies of microwave emission characteristics. 
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